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Abstract. A few tight correlations linking several properties of gamma-ray bursts (GRBs), namely the spectral peak energy, 
the total radiated energy, and the afterglow break time, have been discovered with pie-Swift GRBs. They were used to constrain 
the cosmological parameters, together with type-la supemovae. However, the tightness of these correlations is a challenge to 
ORB models. We explore the effect of adding Swift bursts to the Ghirlanda and Liang-Zhang relations. Although they are both 
still valid, they become somewhat weakened mostly due to the presence of significant outliers, which otherwise are apparently 
normal GRBs so difficult to distinguish. The increased dispersion of the relations makes them less reliable for purposes of 
precision cosmology. 
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1. Introduction 

Gamma-ray bursts (GRBs) release a huge amount of energy 
in a short amount of time, making them the brightest sources 
in the gamma-ray sky. Following these explosions, a residual 
emission, the afterglow, is observed at X-ray and, often, at UV- 
optical-IR wavelengths (see van Paradijs et al. 2000 and Zhang 
& Meszaros 2004 for reviews). Given the involved energies, 
GRBs represent the brightest objects in the Universe and can 
therefore be detectable up to very high redshifts. In contrast to 
type-la supernovae, which are good cosmological rulers up to 
redshift z 1 (e.g. Perlmutter et al. 1999; Riess et al. 2004), it 
has been speculated that GRBs are promising standard candles 
for a high redshift universe, and attempts to use GRB properties 
to indicate distances have been pursued in the past. 

Based on samples of varying size, the isotropic-equivalent 
burst luminosities were found to be correlated with the tem- 
poral variability (smooth bursts are intrinsically less luminous; 
Fenimore & Ramirez-Ruiz 2000; Reichart et al. 2001 ; Guidorzi 
et al. 2005) and anti-correlated with the spectral time lags (the 
times elapsing between the light curve at higher and lower en- 
ergies; NoiTis et al. 2000). There were also encouraging re- 
sults when considering spectral properties. Amati et al. (2002) 
found a correlation between the rest-frame peak energy E'^ — 
(1 + z)Ep and the isotropic-equivalent energy released during 
the prompt phase -Eiso- This correlation, initially based on 10 
GRBs with known redshift, was later updated and confirmed 
using a sample of 41 GRBs (Amati 2006, 2007; see however 
Band & Preece 2005; Nakar & Piran 2005). The Amati rela- 
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tion has been used for a number of statistical studies to de- 
rive the redshift of GRBs with known peak energy and fluence 
(e.g. Atteia 2003; Liang et al. 2004). However, the relatively 
large scatter results in poorly constrained redshifts (especially 
for high values). This indicates that the Amati relation is not 
tight enough for cosmological studies (e.g. Li 2006). To date, 
there are two known outliers of the Amati relation, namely 
GRB 980425 and GRB 031203. These bursts are peculiar in 
several respects, since they are very close (z < 0.1), under- 
luminous with respect to most other events, and they have pe- 
culiar afterglows (e.g. Woosely & Bloom 2006). Driven by the 
observation that XRF060218 showed a strong spectral evolu- 
tion during its prompt emission phase (Campana et al. 2006), 
and that the mean peak energy is consistent with the Amati 
relation (Amati et al. 2007), Ghisellini et al. (2006) speculate 
that these two outliers may also be consistent with the Amati 
relation, if a strong spectral evolution took place. 

A significant step forward was performed by Ghirlanda et 
al. (2004a). Since GRBs are believed to originate inside jets 
(with half-opening angle i9), they considered the collimation- 
corrected energy Ej — (1 — cos i})Eiso instead of the isotropic- 
equivalent value. Remarkably, based on a set of 15 GRBs (now 
expanded to 18; Nava et al. 2006), E^ is much more tightly 
correlated with E^ than is Eiso (see also Friedman & Bloom 
2005). With such tightness, Ghirlanda et al. (2004b) combined 
15 GRBs with data from type-la supernovae and were able to 
constrain the cosmological parameters (see also Xu et al. 2005 
and Schaefer 2007 and references therein). To infer the jet colli- 
mation degree, i) is computed from the time at which the after- 
glow light curve steepens (Rhoads et al. 1999; Sari et al. 1999; 
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Chevalier & Li 2000). There are two possible cases, depend- 
ing on the density profile of the medium surrounding the GRB 
source, which is usually assumed either homogeneous (inter- 
stellar medium, ISM) or wind-shaped (n = Ay. r^"^, wind). In 
particular, we have (Sari et al. 1999; Chevaher & Li 2000) 
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where tb,d is the break time in days, the radiation efficiency, 
-E'iso,52 = i'iso/(10^^ erg), no is the ambient particle density 
in cm^'^, and A* = 3 x 10'^'^ cm^^ is the wind constant. 

The Ghirlanda relation is tight in both cases, with a log- 
arithmic dispersion of 0.10 dex (0.08 dex) in the ISM (wind) 
case; in comparison, the Amati relation has a dispersion of 
0.4 dex. It is not straightforward to compare GRB 980425 and 
GRB 031203 with the GhManda relations, since their collima- 
tion degree is unknown, but it looks likely that these events are 
outliers for them as well. We also note that the Ghirlanda re- 
lation is intimately connected with the jet break interpretation. 
The standard jet break theory predicts a variation in the decay 
index of Aa 1 and a decay index after the break a-i ^ 2. 
These predictions are loosely satisfied by the pie-Swift sample. 
By taking the full analysis of Zeh, Klose & Kann (2006), we 
find that only 3 GRBs of the 13 pre-Swift bursts available (5 
GRBs do not have enough data) have a Aa ^ 1 at 1 <t level. 
Allowing a 3 (T level we have 9 GRBs consistent with this con- 
straint, leaving out 4 GRBs that do not satisfy these relations. In 
addition, 2 bursts definitely have (3 a) a decay which is shal- 
lower than 2. Several different explanations can address this 
issue as the uncertainty of sideways expansion physics (e.g. 
Panaitescu et al. 2005), different electron energy distributions 
(p value) for different GRBs, change in the microphysical pa- 
rameters (e), new energy injections, changing external medium, 
etc. 

The need to estimate the jet angle through Eq. ([T]i some- 
how makes the Ghirlanda relations model-dependent. Recently, 
a lot of discussion has arisen concerning the interpretation of 
breaks in the afterglow light curves of GRBs discovered by 
Swift (e.g. Panaitescu et al. 2006; Sato et al. 2007; Willingale 
et al. 2007; Zhang 2006). For example, Panaitescu et al. (2006) 
present six optical/X-ray afterglows whose light curves clearly 
show a chromatic behavior, which cannot be due to a geomet- 
ric effect. Liang & Zhang (2005) have previously pointed out 
that a purely phenomenological correlation (i.e. without any 
assumption on the nature of the breaks) exists between E'iso, 

^ t^/(l + z), and E'p, in the form Ei^o oc E'^'^fy^^. It should 
be noted that Liang & Zhang (2006) argue that the breaks in 
optical and X-ray bands may have different physical origins. 
They explicitly took th from optical light curves (actually, this 
was also the case for most of the Ghirlanda sample). Nava et al. 
(2006) point out that the Ghirlanda and Liang-Zhang relations 
are mutually consistent, provided that y = —1 (as their data 
support). 

During the past two years, the GRB field has been revo- 
lutionized by Swift. More than 200 long-duration GRBs have 



been discovered and a secure redshift has been obtained for 
more than 50 (up to February 2006). Swift GRBs with known 
redshift are at a higher mean redshift than before (Jakobsson 
et al. 2006). GRBs at larger redshift might have different char- 
acteristics maybe related to metallicity's impact on these cor- 
relations. Obviously it is very relevant to see how these new 
GRBs behave with respect to the GhManda and Liang-Zhang 
relations. To this aim, three parameters are needed for each 
event, namely the redshift, the spectral peak energy, and the 
afterglow break time. The burst alert telescope (BAT) onboard 
Swift, which detects the bursts, covers a limited spectral en- 
ergy range (15-350 keV), so that Ep can only be determined 
in very few cases. However, some GRBs discovered by Swift 
can also be observed by other instruments, in particular by the 
Wind-Konus experiment, which has a useful energy range of 
18-1160/21-1360 keV, or by HETE-2, Suzaku, and RHESSI. 
In this paper, we present a systematic search for all Swift bursts, 
which have the measurement of the three parameters needed for 
comparing them with the GhManda and Liang-Zhang relations. 

When this paper was submitted, we became aware of a sim- 
ilar paper (Ghirlanda et al. 2007). We have not carried out a de- 
tailed comparison between the two works but highlight a few 
differences in the following sections. 



2. Sample selection 

We collected all long-duration GRBs detected by Swift that 
have a secure (spectroscopic) redshift z and a published peak 
energy Ep. We found 19 events satisfying this criterium. We 
thus inspected their afterglow light curve looking for breaks. 
Due to the unclear multi-wavelength behavior of these breaks, 
we did not look at X-ray or radio data. This is mainly dic- 
tated by the fact that, when the Ghirlanda relation (and simi- 
lar works) was found, no simultaneous optical and X-ray light 
curve existed. This choice (sticking to optical light curves) also 
guarantees an unbiased extension of the GhManda relation. In 
any case, when describing the GRBs below, we outline when 
differences from the standard jet model are occurring. We end 
up with eight GRBs that have a reliable measurement of ib 
(GRB 060614 is not used however). For eight more, we can 
only provide a lower limit to this parameter For the other three 
there is not published information (yet) to constrain it. To com- 
pute the beaming-corrected energy, we first need the isotropic- 
equivalent value. Recently, Amati (2006, 2007) listed the val- 
ues of Eiso and E'^ for 18 Swift bursts, and we adopted his val- 
ues whenever available. In other cases, we resorted to values 
published in the GCN^] (see below), computing a bolometric 
correction when necessary. Following Ghirlanda et al. (2004), 
for all bursts we assumed an external density no = 3 cm^'^ 
and a wind parameter A^, — 1 (for both quantities we also in- 
troduced an error of Alogj^g n{A^) = 0.5). When dealing with 
lower limits in the break times, we conservatively assumed a 
50% error on that limit. In addition, we consider the errors in 
logarithmic space by taking the left/right error closer to the best 
fit line (this has been done also for GRBs in the Nava's sam- 
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Table 1. Sample of Swift GRBs with known z, Ep, and fb- 



GRB 


z 


E' 
p 




log -Biso 






log M"' 


log M"' 


Pred. t^,*"' 


Pred. t^^^ 






(keV) 


(d) 


(erg) 


(deg) 


(deg) 


(erg) 


(erg) 


(d) 


(d) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


050525A 


0.606 


127±10 


0.3±0.1 


52.53±0.03 


4.0±0.7 


3.7±1.0 


49.91±0.15 


49.86±0.16 


5+"-i 


5+"1 

" 0.2 


050820A 


2.612 


1325±270 


15±8 


53.92±0.09 


8.5±2.0 


3.7±1.1 


51.96±0.17 


51.23±0.16 


4.2+?-? 


— A.o 


050922C 


2.198 


415±111 


0.11 ±0.03 


52.79±0.17 


2.0±0.3 


2.1±0.6 


49.55±0.23 


49.62±0.27 


*^ — 2.5 


^ —3.2 


05 11 09 A 


2.346 


539±381 


0.60±0.10 


52.70±0.22 


3.7±0.4 


3.4±0.6 


50.02±0.26 


49.94±0.32 


ii+f 

— y 


12+f 

— 9 


060124 


2.297 


636±162 


1.13±0.09 


53.62±0.06 


3.6±0.5 


2.3±0.6 


50.93±0.09 


50.54±0.15 


— 0.8 


2.0+?-^ 

— l.U 


060206 


4.045 


381±98 


0.6±0.2 


52.76±0.08 


3.1±0.6 


2.9±0.8 


49.93±0.17 


49.88±0.19 


7.6+|-^ 


8.o+!-5 


060526 


3.21 


105±21 


2.77±0.30 


52.41±0.05 


5.5±0.4 


3.9±0.5 


50.23±0.08 


50.07±0.11 


0.3+°-? 

— U.l 


0.3+°? 

— 0.1 


060614 


0.125 


55±45 


1.39±0.04 


51.40±0.22 


11.1±0.9 


11.5±1.8 


49.67±0.25 


49.70±0.31 


'-'•-^—0.2 


" 0.2 


050318 


1.44 


115±25 


>0.26 


52.41±0.03 


>3.3 (0.7) 


>3.5 (0.6) 


>49.63 (0.22) 


>49.68 (0.20) 


0.8+r, 

^^—0.3 


0.9+"« 

" — 0.4 


050401 


2.90 


501±53 


>13 


53.61±0.08 


>8.5 (2.0) 


>4.1 (1.2) 


>51.66 (0.16) 


>51.03 (0.15) 


1.4+°t 


1.5+07 


0504 16A 


0.653 


25.1±4.2 


>4.7 


51.08±0.07 


>16.7 (3.8) 


>17.0 (5.0) 


>49.70 (0.15) 


>49.72 (0.15) 


6+°" 


6+°'® 


050603 


2.821 


1333±107 


>2.5 


53.85±0.03 


>4.3(1.0) 


>2.4 (0.7) 


>51.30 (0.22) 


>50.79 (0.20) 


5 3+1-^ 




060418 


1.489 


572±114 


>5 


52.95±0.05 


>8.5 (2.0) 


>5.3 (1.6) 


>5 1.00 (0.23) 


>50.59 (0.20) 


5.31?:^ 


5.61^2:? 


060927 


5.6 


473±116 


>0.16 


52.92±0.09 


>1.6 (0.4) 


> 1.8 (0.5) 


>49.54 (0.25) 


>49.61 (0.23) 


loti 


11+12 


061007 


1.261 


902±43 


>1.74 


54.00±0.04 


>5.3(1.2) 


>2.6 (0.8) 


>51.47(0.15) 


>50.90 (0.14) 


1 o+o-^' 

^•"-0.2 


1 1+0.8 
^•^-0.5 


061121 


1.314 


1288±153 


>3.5 


53.36±0.04 


>6.8(1.6) 


>3.9(1.1) 


>51.21 (0.23) 


>50.73 (0.20) 


t s+o-'^ 

i-°-0.5 


1 Q+1-5 



(1) - GRB name. (2) - redshift. (3) - Rest-frame spectrum peak energy. (4) - Observed break time in the optical band and in the observer 
frame. (5) - Isotropic-equivalent bolometric energy. (6), (7) - Jet opening angle computed for homogeneous and wind environments. (8), (9) - 
Beaming-corrected energy for the homogeneous and wind cases. (10), (1 1) - Predicted break time assuming the Ghirlanda relation as computed 
by Nava et al. (2006) after correcting for a small typo in the published relation (Nava, private communication). 
GRB060614 is not used in the analysis described in the text, but it is reported here for completeness. 



pie). Below we summarize their properties and all the values 
are reported in Table [T] 

• GRB 050318 was observed by Swift and the BAT was 
able to determine the peak energy (Perri et al. 2005). From the 
UVOT data, one can determine a lower limit on the time of a 
possible break in the optical light curve of 0.26 d (Still et al. 
2005). In the X-ray, a break around 0.2 d is also seen (Perri et 
al. 2005). 

• GRB 050401 was observed both by 5w//f-BAT and Wnd- 
Konus. We take for the average of the two values presented 
by Golenetskii et al. (2005a). In the optical, the well-monitored 
light curve (Rykoff et al. 2005; Watson et al. 2006) showed a 
slow, uninterrupted decay up to 13 d after the burst. In the X- 
ray band, the afterglow was observed up to 11.6 d after the 
trigger (De Pasquale et al. 2005), revealing a break at 0.06 ± 
0.01 d (a value inconsistent with the Ghirlanda relation; De 
Pasquale et al. 2005). Ghirlanda et al. (2007) indicate a possible 
break in the optical light curve at 1.5 ± 0.5 d. We note that 
a simple power-law fit is able to reproduce all the observed 
optical data (Xrod ~ ^-^ ^^'-^ ^8 degrees of freedom, with a 
null hypothesis probability of 17%). We also tried to fit the 
same data with a smoothly-joined power-law without finding 
any statistical evidence of a break. 

• XRF050416A has a well-defined X-ray Hght curve 
(Mangano et al. 2007a), showing the common steep-flat-steep 
behavior (Nousek et al. 2006; O'Brien et al. 2006). No addi- 
tional breaks are visible between 0.017 and ^ 42 d after the 
GRB (Sakamoto et al. 2006; Mangano et al. 2007a). Thanks 
to its softness it was possible to compute the peak energy with 
BAT (Sakamoto et al. 2006). In the optical there are also no 
signs of a break in the light curve up to at least ~ 4.7 d (see the 



/-band light curve presented by Holland et al. 2007; see also 
Soderberg et al. 2007). Ghirlanda et al. (2007) fitted the optical 
light curve, including a XRF 0602 1 8 -like light curve (Campana 
et al. 2006). In their fit they find a break at 1.0 ± 0.7 d. We note 
that the SN presence is suggested by Soderberg et al. (2006) 
and they do not require the presence of a break in the light 
curve in their fit with the SN modeled with a SN1998bw tem- 
plate. 

• GRB 050525A was a bright burst observed by Swift and 
W/nt/-Konus (Blustin et al. 2005; Golenetskii et al. 2005b). It 
has already been included in the sample of Nava et al. (2006). 
Several groups (Blustin et al. 2005; Mirabal et al. 2005; Delia 
Valle et al. 2006a) report a break in the optical about 0.3 d after 
the GRB. We adopt here — 0.3 ± 0.1 d, which is the average 
of the values found by these authors (see also Nava et al. 2006). 
In the X-ray band, a break was also detected, at about the same 
time as the optical one, suggesting a jet origin for this afterglow 
(Blustin et al. 2005). Sato et al. (2007), however, question this 
interpretation, since the post-break decay indices were different 
in the optical and X-ray bands, contrary to model expectations, 
and the decay index after the break is too shallow in the X-ray 
and optical bands. 

• The X-ray light curve of the GRB 050603 afterglow was 
a simple power law between 0.4 and 7 d after the GRB (Grupe 
et al. 2006). The peak energy was estimated through Wind- 
Konus observations (Golenetskii et al. 2005c). UVOT observed 
the optical afterglow as well, showing no breaks within the in- 
terval 0.4-2.5 d (Grupe et al. 2006). A slope a = 1.8 ± 0.2 
was found by UVOT, consistent with the XRT one. Given the 
relatively steep value, it is unclear whether a break should have 
appeared before or after the Swift observations. The Ghirlanda 
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relation would predict a jet break time around 4-7 d, so we 
consider th > 2.5 d. 

• GRB 050820A was a bright burst with large Ei^o and 
(Wind-Konus, Cenko et al. 2006). The X-ray light curve of 
GRB 050820A shows the usual triple power-law behavior. A 
lower limit to any further break in the X-rays can be set to 

~ 17 d. The optical light curve tracks the X-ray one, but 
thanks to late-time HST images a break is detected at very late 
times, and its timing is difficult to constrain. Using the data 
from Cenko et al. (2006) and allowing for a late-time temporal 
slope in the range 1.5 < a < 3, one gets 7 d < < 23 d. This 
is consistent with the value ib = 18 ± 2 d found by Cenko et 
al. (2006) which imposed a = 2.34 ± 0.06. 

• The X-ray and optical light curves of GRB 050922C have 
been presented by Panaitescu et al. (2006). Both the X-ray 
and optical light curves show a break, but at different times. 
A multi-band fit of the UVOT optical data provides a break 
time at 2.7 ± 0.7 hr (Li et al. 2005; see also Andreev et al. 
2005). The currently published data show no further breaks up 
tot ~ 1.2 d. The peak energy was estimated by HETE-2 ob- 
servations (Crew et al. 2005), even if Wind-Komxs has failed to 
reveal an exponential cut-off (Golenetskii et al. 2005d). 

• GRB 05 11 09 A was discovered by Swift and also observed 
by Wind-Konus (Golenetskii et al. 2005e). The sparse optical 
data constrain the break time to be <b > 0.64 d (Pavlenko et al. 

2005) . More recently Yost et al. (2007), suggested the presence 
of a break occurring at 0.6 ± 0.1 d in the optical light curve. 
This break however does not coincide with an X-ray break and 
is interpreted as a cooUng break. 

• GRB 060124 is an astonishing burst. Swift was triggered 
by the precursor, allowing to study in detail the prompt phase 
with CCD X-ray spectroscopy for the first time (Romano et 
al. 2006). In the X-ray band, a break in the light curve occurs 
at 1.21 ± 0.17 d. At variance with many other Swift GRBs, a 
break in the optical light curve was observed nearly simulta- 
neously with the X-ray break at 1.13 ± 0.09 d (D. A. Kann, 
priv. comm.; see also Curran et al. 2007). The spectral peak en- 
ergy was computed by simultaneously fitting Wind-Konus and 
BAT data (Romano et al. 2006); note that the value reported by 
HETE-2 (Lamb et al. 2006) is slightly higher, but still consis- 
tent with the adopted value. Also in this case, the X-ray decay 
index after the break is too shallow for a jet break. 

• The late optical and X-ray light curves of GRB 060206 
have been published by Monfardini et al. (2006) and Stanek 
et al. (2006). In both bands the light curve is complex, but a 
clear break in the optical curve can be detected at fb = 0.6 d. 
The light curve then goes on uninterrupted until at least 2.3 d 
(Stanek et al. 2006; Monfardini et al. 2006). We take the Ep 
and Eiso from the value measured by Swift-BKT (Palmer et al. 

2006) . 

• The X-ray Ught curve of GRB 060418 shows several 
bright flares. On the other hand, the optical/NIR light curve is 
smooth (Molinari et al. 2007), with no breaks between ~ 150 s 
and 5 d. The Ep and Eiso values have been taken from the 
Konus data (Golenetskii et al. 2006a). 

• GRB 060526 was discovered by Swift, showing two 
strong flares just after the main event. A simultaneous break in 
the optical and X-ray light curves is observed at ^ 2.8 d (Dai et 



al. 2007). The bolometric source flux and the peak energy were 
taken by Schaefer (2007). 

• GRB 060614 was a remarkable low-redshift burst that 
did not show supernova signatures down to very deep hmits 
(Delia VaUe et al. 2006b; Fynbo et al. 2006a; Gal- Yam et al. 
2006). The early-time light curve was complex, but, remark- 
ably, VLT and Swift-XRT data show an achromatic break at 

1.3 d (Mangano et al. 2007b). The peak energy of this 
event is constrained by the Wind-Konus data (Golenetskii et 
al. 2006b), which present a time-resolved analysis. Amati et al. 
(2007) have estimated the average Ep to he in the range 10- 
100 keV range. Since the BAT spectrum shows no deviation 
from a power-law (Gehrels et al. 2006) in the range 15-150 
keV, we further assumed Ep < 50 keV. We took the value of 
Eiso from Amati et al. (2007). Given the lack of a bright super- 
nova, however, there is an ongoing debate about whether this 
GRB is related to a massive progenitor or to a merging of two 
compact objects (Gehrels et al. 2006; Zhang et al. 2006). Since 
the Amati and Ghirlanda relations hold for long GRBs, we do 
not consider this burst in the following. 

• GRB 060927 was a high-redshift burst at z = 5.6 (Fynbo 
et al. 2006b). The Swift observations allow Ep to be con- 
strained inside the BAT range. We computed the rest-frame 1- 
10000 keV fluence by employing the spectral parameters pro- 
vided by Stamatikos et al. (2006) and found Ei^o = (8.4 ± 
1.6) X 10^^ erg. The optical data only imply a Umit t^, > 0.1 d 
(e.g. Antoniuk et al. 2006). The X-ray light curve shows a 
break at t = 0.05 d (Troja et al. 2006). 

• The X-ray and optical light curves of GRB 061007 con- 
sist of a single power law starting soon after the burst. Current 
limits show tb > H and 1.7 d after the GRB in the X-ray 
and optical bands, respectively (Mundell et al. 2007; Schady 
et al. 2007). Schady et al. (2007) suggest that a break might 
have occurred before the start of the XRT and UVOT obser- 
vations (72 s after the burst). However, ROTSE observations 
starting 26 s after the burst trigger imply a rising flux between 
the two measurements, making this hypothesis unlikely. Thus 
we chose tb > 1-7 d. We took Ep and ii^iso from Golenetskii 
et al. (2006c), whose results are consistent with the RHESSI 
values (Wigger et al. 2006). We note that Suzaku-WAM found 
a larger Ep (Yamaoka et al. 2006). 

• GRB 061 121 showed a break at X-ray wavelengths at 
3 ks after flie GRB (Page et al. 2006). In the optical band, there 
is a possible rebrightening at ^ 3 ks, and the decay is then 
regular up to 3d (Halpern & Armstrong 2006; Efimov et 
al. 2006). Determining the location of a break is complicated 
by the presence of a bright host galaxy (Malesani et al. 2006; 
Cobb 2006). We took the average of the Ep values as found by 
Wind-Konus (Golenetskii et al. 2006d) and RHESSI (BeUm et 
al. 2006). 

• There are three more GRBs with measured z and Ep, 
namely GRB 0601 15, GRB 060218, and GRB 060707. For 
these bursts, it is however impossible to measure a break time, 
due to the small amount of available data or to the presence 
of SN2006aj in the case of GRB 060218 (e.g. Campana et al. 
2006; Pian et al. 2006). We encourage observers with avail- 
able data to provide a measurement for the break time of these 
bursts. 
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The values that we adopted for the optical break times are 
strictly compliant to the selection made by Ghirlanda et al. 
(2004a). The assumed error budget is conservative. By taking 
the limits from the Swift X-ray light curves (which often ex- 
tend much further in time), we would find even tighter limits, 
as recently shown by Sato et al. (2007) for a subsample of our 
GRBs. 



3. Analysis and results 

We added the new Swift bursts to the sample of 18 events pre- 
sented by Nava et al. (2006). We note that GRB 050525A is 
common to both samples. There are eight Swift bursts with a 
measurement of i^iso, Ep and th (GRB 0606 14 is not used how- 
ever), plus eight events with a lower limit on tb- In Fig. [T]and|2] 
we show the Swift bursts in the E'^ vs. E-^so plane. 

We repeated the analysis caiTied out by Nava et al. (2006) 
on their sample using the routine ^fexy of Press et al. (2003) 
and fit the Ghirlanda relation to their data for both the ISM 
and wind cases. We obtained consistent results (see Table 2), 
adopting the value of no = 3 cm^^ for the circumburst density 
and Tj-f — 20% for the gamma-ray efficiency for all events. 

Before starting the statistical analysis, we observed that 
five bursts with lower limits on the break time lie at the left 
of the best-fit Ghirlanda relation (GRB 050318, GRB 050603, 
GRB 060418, GRB 060927 and GRB 061121, see Figs. [Hand 
|2|i. Since they have only lower limits on the corrected energy 
and are consistent with the relation, they are not considered in 
the fit. 

We were left with 27 GRBs. As can be seen in Figs. 1 and 
2, there are significant outliers to the relation. Adopting the 
routine ^fexy, we derived a reduced = 3.82 (25 dof) for 
the homogeneous ISM relation (null hypothesis probability of 
4 X 10-1°) and a reduced = 3.09 (25 dof) for a wind ISM 
relation (null hypothesis probability of 3 x 10"^). We note that 
our limits are conservative, since we had treated the lower lim- 
its on as actual measurements. If the breaks occurred sig- 
nificantly later than these limits, many of the bursts marked by 
arrows in Figs. [T]and|2] would move rightward in the plot, even 
more increasing the scatter in the coiTelation. 

Given the concerns about using the Ghirlanda relation 
in terms of the jet interpretation, we also fitted the model- 
independent Liang-Zhang relation to the bursts in the sam- 
ple. Again, by including only the events listed in the Nava et 
al. (2006) sample, we found consistent results with what they 
found, albeit with a slightly larger (see Table 2). Using the 
data as listed in the Tables 1 and 3 of Firmani et al. (2006), 
we find a significantly larger x^, which is entirely due to 
GRB 050525A (variation in the reduced from 0.83 to 2.46). 
We then added the Swift bursts to the sample, and again found 
that the relation worsened significantly. Moreover, the expo- 
nent y in the Liang-Zhang relation is now y — —0.55 ± 0.10 
(68% confidence level), significantly different from the value 
y = —1 needed to make the Ghirlanda and Liang-Zhang rela- 
tions mutually consistent (Nava et al. 2006). 



4. Discussion and conclusions 

We have investigated in detail the addition of Swift bursts to the 
Ghirlanda and Liang-Zhang relations. Given the narrow energy 
band of the BAT instrument, only a small fraction of the > 200 
GRBs so far discovered by Swift precisely determine the peak 
energy needed for the comparison with the Ghirlanda correla- 
tion. Putting together data from different instruments (mostly 
W/«(i-Konus, and occasionally HETE-2, RHESSl, and Suzaku) 
and including the few bursts for which a prompt measurement 
of Ep was provided by BAT, we find a total of 19 GRBs. Of 
these, eight have a determination of the optical light curve 
break (seven are used), eight have a lower limit (three are used), 
and three not have enough (published) optical data to constrain 
it. 

With the inclusion of these nine new GRBs (one was al- 
ready in the Nava sample), the Ghirlanda correlation suffers 
from a worsening. Fitting the data, the reduced x^ goes from 
1.1 (1.1) to 3.8 (3.1) in the case of a homogeneous (wind) 
medium (see Table 2). The correlation is clearly still present, 
but it is not as tight as suggested in previous works. Its use for 
cosmological purposes is thus weakened. Actually, Spearman's 
and Kendall's (tau) correlation tests give a relatively weak 
probability of correlation on the Swift sample (seven objects). 
With these tests we obtain correlation probabilities of ~ 20% 
for both homogeneous and wind cases, to be compared with 
10"'' - 10"* for the pre-Swift sample. 

In Table 1 we also report the predicted break times based on 
the Ghirlanda relation in order to see if there are correspond- 
ing breaks in the X-ray light curves. On average the predicted 
breaks do not occur in concordance with X-ray breaks. This 
means that even when considering breaks in the X-ray light 
curves, we are not able to reconcile the observed times with the 
ones predicted by the Ghirlanda relation. 

The optical light curves of bursts discovered by Swift 
are significantly more complex than previously thought (e.g. 
Stanek et al. 2006; Gal-Yam et al. 2006). One might there- 
fore speculate that the breaks observed in some of the strong 
outliers are different from those responsible for the existence 
of the Ghirlanda correlation. For example, GRB 050922C and 
GRB 060206 had breaks both in the optical and in the X- 
ray light curves, but at different times (Panaitescu et al. 2006; 
Stanek et al. 2006). The use of these breaks (so-called chro- 
matic breaks) for the Ghirlanda relation is not appropriate in 
principle, since a jet break produces an achromatic steepening. 
Only GRB 060124, GRB 050525A and GRB 060614 show ev- 
idence for an achromatic break at a known redshift. Recently, 
several papers suggest a different origin of the X-ray and op- 
tical emission (Uhm & Belobodorov 2007; Grenet, Daigne & 
Mochkovitch 2007; Ghisellini et al. 2007), so that the need for 
achromatic breaks might be relaxed. 

Optical breaks not coincident with X-ray breaks might 
have been missed before the launch of Swift due to the lack 
of X-ray coverage. For the outliers lying above (below) the 
correlation, a further break at later (earlier) times (the true "jet- 
break") would be necessary in order to reconcile them with 
the Ghirlanda relation. This does not always look like a vi- 
able solution. For example, considering outliers in the upper 
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Table 2. Fit results 
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Correlation 


Nava [our] analysis 


Swift data 


Swift data achromatic breaks 


Swift data pure breaks 




X?cd (dof) 


xLd (dof) 


X?cd (dof) 


X?od (dof) 


Ghirlanda ISM 


1.40 [1.42] (16) 


3.82 (25) 


2.35 (22) 


2.33 (20) 


Ghirlanda wind 


1.13 [1.13] (16) 


3.09 (25) 


2.00 (22) 


1.72 (20) 


Liang & Zhang* 


1.49 [1.81] (15) 


4.56 (24) 


2.44 (21) 


2.37 (19) 



* The logarithmic errors in this case have been evaluated as the mean of the lower and higher errors. 



part of the E'^-E-y plane, GRB 060206 showed a steep decay 
(a ^ 2) after 0.6 d (IVIonfardini et al. 2006; Stanek et al. 
2006), similar to many of the breaks used to build the original 
Ghirlanda relation. Even if no observations are available at late 
times, a further break appears unlikely. Late observations were 
probably secured (but not published) for GRB 050922C, which 
would be particularly interesting in this respect. Concerning the 
bursts lying on the right side of the correlation, early observa- 
tions are available for GRB 050401 (Rykoff et al. 2005; Watson 
et al. 2006), which shows a remarkable, unbroken decay start- 
ing 35 s after the burst. Thus, at least for this burst, we have 
strong "evidence for absence" of a break that might have rec- 
onciled GRB 050401 with the Ghirlanda relation. 

In addition, considering only bursts with achromatic 
breaks and upper limits (i.e. we disregard GRB050922C, 
GRB051109A and GRB 060206; we retain GRB050820A 
since we do not have information from the X-rays), the re- 
sults of the fits are improved, with a reduced of 2.35 in the 
case of a homogeneous medium (null hypothesis probability of 
4 X 10^^) and with a reduced — 2.00 (null hypothesis prob- 
ability of 0.003) in the case of a wind medium. As described 
above, concerns have also been expressed on GRB 050525A 
and GRB 060124 since their post break decay index is too flat 
and the decay variation is less than one. We thus also consider 
the case of only GRB 050820A and GRB 060526 and the upper 
limits (pure breaks in Table 2). In this case we obtain a reduced 

of 2.33 for a homogeneous medium (null hypothesis proba- 
bility of 7 X 10^*) and — 1.72 (null hypothesis probability 
of 0.02) for a wind medium. 

One might speculate on whether the Swift bursts have pe- 
culiar properties. In this respect, it is interesting to note that 
neither GRB 050922C nor GRB 050603 (the two strongest out- 
liers of the relation) have any particularity, in terms of their 
duration, spectral properties, and fluence in the Swift sample. 
A distinguishing property of Swift bursts as a sample is their 
average measured higher redshift (e.g. Jakobsson et al. 2006). 
The average redshift of the Swift burst introduced in our sam- 
ple is (z) — 2.2, to be compared with (z) — 1.4 for the events 
in the sample of Nava et al. (2006). This might hint at a pos- 
sible evolutionary effect, possibly due to the the properties of 
the jet propagation inside the progenitor stars and then in their 
environment during the afterglow phase. The current sample is, 
however, still too limited to draw any strong conclusion. 
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49 50 51 52 
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Fig. 1. The rest-frame spectral peak energy E'p and the colHmation coiTected energy E^. This energy has been calculated assum- 
ing a homogeneous density medium with a value of rto = 3 cm^"^. The open circles mark the pTS-Swift bursts used to derive the 
Ghirlanda relation in Nava et al. (2006). Filled circles represent Swift GRBs with known peak energy, redshift, and break time in 
the optical. The arrows mark Swift bursts with only lower limits in the break time tb- Those events marked with stars he at the left 
of the best-fit correlation and are consistent with it, so they have not been included in the fit. Events with a big circle around are 
the two GRBs with the optical break not coincident with the X-ray break, which have been discarded in the additional analysis. 
The point marked with a diamond is GRB 060614, shown for comparison purposes. The dashed line shows the best-fit power law 
model for the events in the Nava et al. (2006) sample, obtained by accounting for the errors on both coordinates. The solid line 
shows the same fit to the previous sample with the addition of all Swift bursts (except those marked with stars). 
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Fig. 2. As in the previous figure but for a wind medium with density parameter A* = 1. 



